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Abstract 



This article deals with the flow of Newtonian fluids through axially-symmetric 
corrugated tubes. An analytical method to derive the relation between volumetric 
flow rate and pressure drop in laminar flow regimes is presented and applied to a 
number of simple tube geometries of converging-diverging nature. The method is 
general in terms of fluid and tube shape within the previous restrictions. Moreover, 
it can be used as a basis for numerical integration where analytical relations cannot 
be obtained due to mathematical difficulties. 



1 INTRODUCTION 
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1 Introduction 

Modeling the flow through corrugated tubes of various geometries is an important 
subject and has many real-life applications. Moreover, it is required for modeling 
viscoelasticity, yield-stress and the flow of Newtonian and non-Newtonian fluids 
through porous media [1-3]. There are many previous attempts to model the flow 
through capillaries of different geometries. However, they either apply to tubes of 
regular cross sections [4, 5] or deal with very special cases. Most of these studies 
use numerical mesh techniques such as finite difference and spectral methods to 
obtain numerical results. Illuminating examples of these investigations are Kozicki 
et al. [6], Miller [7], Oka [8], Williams and Javadpour [9], Phan-Thien et al. [10, 11], 
Lahbabi and Chang [12], Burdette et al. [13], Pilitsis et al. [14, 15], James et al. 
[16], Talwar and Khomami [17], Koshiba et al. [18], Masuleh and Phillips [19], and 
Davidson et al. [20]. 

In the current paper we present an analytical method for deriving the relation- 
ship between pressure drop and volumetric flow rate in corrugated tubes of circular 
but varying cross section, such as those depicted schematically in Figure 1. We also 
present several examples of the use of this method to derive equations for New- 
tonian flow although the method is general and can be applied to non-Newtonian 
flow as well. In the following derivations we assume a laminar flow of a purely- 
viscous incompressible fluid where the tube corrugation is smooth and limited in 
magnitude to avoid problematic flow phenomena such as vortices. 



2 CORRELATION BETWEEN PRESSURE DROP AND FLOW RATE 






Figure 1: Profiles of converging-diverging axisymmetric capillaries. 

2 Correlation Between Pressure Drop and Flow 
Rate 

A Newtonian fluid with a viscosity fj, satisfy the following relation between stress 
r and strain rate 7 

T = fi'y (1) 

The Hagen-Poiseuille equation for a Newtonian fluid passing through a cylin- 
drical pipe of constant circular cross section, which can be derived directly from 
Equation 1, states that 

Q = -^ (2) 

where Q is the volumetric flow rate, r is the tube radius, P is the pressure drop 
across the tube, /i is the fluid viscosity and x is the tube length. A derivation of 
this relation can be found, for example, in [1, 21]. On solving this equation for P, 



2.1 Conical Tube 
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the following expression for pressure drop as a function of flow rate is obtained 

_ 8Qfix 

For an infinitesimal length, 6x, of a capillary, the infinitesimal pressure drop for 
a given flow rate Q is given by 

6P = ^ (4) 

For an incompressible fluid, the volumetric flow rate across an arbitrary cross 
section of the capillary is constant. Therefore, the total pressure drop across a 
capillary of length L with circular cross section of varying radius, r{x), is given by 

p ^ 8Q/i dx 
7^ Jo 

This relation will be used in the following sections to derive relations between 
pressure drop and volumetric flow rate for a number of geometries of axisymmetric 
capillaries of varying cross section with converging- diverging feature. The method 
can be equally applied to other geometries of different nature. 

2.1 Conical Tube 

For a corrugated tube of conical shape, depicted in Figure 2, the radius r as a 
function of the axial coordinate x in the designated frame is given by 

r{x) = a + b\x\ -L/2<x<L/2 (6) 

where 



a = Rmin and b = (7) 



2.2 Parabolic Tube 

Hence, Equation 5 becomes 
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Figure 2: Schematic representation of the radius of a conically shaped converging- 
diverging capillary as a function of the distance along the tube axis. 



2.2 Parabolic Tube 

For a tube of parabolic profile, depicted in Figure 3, the radius is given by 



r{x) = a + bx'^ - L/2 < x < L/2 (12) 



where 



2.2 Parabolic Tube 



9 



r 



R 



-U2 







U2 



Figure 3: Schematic representation of the radius of a converging- diverging capillary 
with a parabolic profile as a function of the distance along the tube axis. 



Therefore, Equation 5 becomes 



and ^ ~ L (-^"^"^ ~ Rmm) 



(13) 



Qfi /"^/^ dx 



P , 

J-L/2 {a + bx'^y 

On performing this integration, the following relation is obtained 



(14) 



P 



8Qfi 



TT 



X 



+ 



5x 



+ 



5x 



5 arctan ( x 



6a{a + bx'^)^ 24a^(a + 6x^)^ 16a^{a + bx'^ 



+ 



L/2 



. -L/2 

(15) 



L 5L 5L 

f-\ A or , 7/-r /r-t\ono ^„ or , 1 / -r / \ 



10 arctan (f 



6a[a + 6(L/2)2]3 24a2[a + 6(L/2)2]2 16a3[a + 6(L/2)2] 1607/2^6 

(16) 



that is 



2.3 Hyperbolic Tube 
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ALQfi 



1 5 5 

+ TTZ^ + 77^^ ^ + 



5arctan ^" 



"^RminRmax ^'^■Rmin^'max ^R^in^max SR^f^y/ Rmax ~ R-n 



(17) 



2.3 Hyperbolic Tube 

For a tube of hyperbolic profile, similar to the profile in Figure 3, the radius is 
given by 

r{x) = Va + hx^ - L/2 < x < L/2 a, 6 > (18) 

where 



and h=[-] (i?L. - RLn) (19) 



Therefore, Equation 5 becomes 



8Q/i /"^/^ dx 



J-L/2{a + bx^y 



IT 



X 



+ 



a.rctan{x^yb/a) 



2a{a + 6x2) 2aVab 



L/2 



-L/2 



(20) 
(21) 



that is 



2.4 Hyperbolic Cosine Tube 
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2.4 Hyperbolic Cosine Tube 

For a tube of hyperbolic cosine profile, similar to the profile in Figure 3, the radius 
is given by 



r(x) = acosh(6x) — L/2 < x < L/2 



(23) 



where 



n 17" 1 / Rmax 

a = Kmin and o = — arccosh 



Hence, Equation 5 becomes 



2 

T 



(24) 



P 



8Qfi 



L/2 



dx 



^ J -L/2 [acosh(6a;)] 



TT 



tanh(6x) [sech^(6x) + 2] 



L/2 
-L/2 



(25) 
(26) 



that is 



P 



37ri?i,„ 



tanh (arccosh ( ) ) sech^ (arccosh ( %^ 



arccosh ( 



2.5 Sinusoidal Tube 



(27) 



For a tube of sinusoidal profile, depicted in Figure 4, where the tube length L spans 
one complete wavelength, the radius is given by 



r{x) = a — 6 cos (kx) 



L/2 < X < L/2 a>6>0 



(2^ 



2.5 Sinusoidal Tube 
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Figure 4: Schematic representation of the radius of a converging- diverging capillary 
with a sinusoidal profile as a function of the distance along the tube axis. 



where 



Hence, Equation 5 becomes 



P 



8Qfi 



L/2 



dx 



^ J -L/2 [a — hcos{kx)] 
On performing this integration, the following relation is obtained 



(30) 



P 



[I] 



(31) 



where 



(6^3 + 9A) ^ /'(A-l)tan('="^ 



arctan 



^11^2 + 4) sin(A;x) 



3(^2 - 1)7/2 

bA sin(A;x) 



- 1 j 6(^2 - l)3[A + cos(A;x)] 
sin(fca;) 



6(^2 - If [A + cos(A;x)]2 3(^2 -l)[A + cos{kx)Y 



(32) 



A 



Rmin Rmax 



(33) 



2.5 Sinusoidal Tube 
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On taking lim _l + I and lim ^- / the following expression is obtained 



Tib^k 



{QA^ + 9A) TT {QA^ + 9A) tt 



3(^2 - 1)7/2 2 3(^2 - 1)7/2 2 



8Qfi{6A^ + 9A) 
'364A;(A2- 1)7/2 



(34) 
(35) 



Since A < —1, P > as it should be. On substituting for A, b and k in the last 
expression we obtain 



Rmax ^min J \ Rmax Flmin 



(36) 



that is 



max \ ^ ^mm max ~l~ Rmin) i^Rmax Rn 

2 TT ( Rmax Rmin ')'^ ^ 



(37) 



It is noteworthy that all these relations (i.e. Equations 11, 17, 22, 27 and 37), 
are dimensionally consistent. Moreover, they have been extensively tested and 
verified by numerical integration. 



3 CONCLUSIONS 
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3 Conclusions 



The current paper proposes an analytical method for deriving mathematical re- 
lations between pressure drop and volumetric flow rate in axially symmetric cor- 
rugated tubes of varying cross section. This method is applied on a number of 
converging-diverging capillary geometries in the context of Newtonian flow. The 
method can be equally applied to some cases of non-Newtonian flow within the 
given restrictions. It can also be extended to include other regular but non-axially- 
symmetric geometries. The method can be used as a basis for numerical integration 
when analytical solutions are out of reach due to mathematical complexities. 



Nomenclature 
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T 



fluid viscosity (Pa.s) 
stress (Pa) 



P 



L 



tube length (m) 
pressure drop (Pa) 



Q 




r 



■max 



tube radius (m) 

maximum radius of corrugated tube (m) 



R. 



■mm 



minimum radius of corrugated tube (m) 



X 
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